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Abstract. - On the basis of a macroscopic ground state population it was argued recently 
that Bose-Einstein condensation should occur in a one-dimensional harmonic potential. We 
examine this situation by drawing analogies to Bosons in a two-dimensional box, where the 
thermodynamic limit is well-defined. We show that in both systems although the ground state 
populations show sharp onsets at the critical temperature, the behaviour of the specific heat is 
analytic, which proves the absence of a phase transition in these systems. 



The experimental study of ultracold trapped Bose gases g, has revived the interest in 
Bose-Einstein condensation in the regime of weak or even vanishing interaction. In these 
experiments the atoms are confined by external forces, which might be modelled by a three- 
dimensional harmonic potential. Recently it was suggested on the basis of results for the ground 
state population that in an effectively one-dimensional potential Bose-Einstein condensation 
should occur Q . In the present paper we place the discussion of non- interacting Bosons in a 
one-dimensional harmonic potential into a larger framework by drawing analogies to particles 
in a two-dimensional box with infinite walls, a system which has, in contrast to the former 
one, the advantage that a well-defined thermodynamic limit exists. Within this framework we 
address the conceptual problem whether a macroscopic ground state population is a sufficient 
indicator for Bose-Einstein condensation to appear or if in contrast, thermodynamic quantities 
have to be considered. In addition to the ground state population we study in particular the 
specific heat, calculated with a continuous density of states, which is permissible for large 
particle numbers. The result is confirmed by an analysis based on the discrete level structure. 
Bose-Einstein condensation Q is usually described within the grandcanonical ensemble 
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where the average number of particles is given by 



N= / AE P {E) — . (1) 

J exp(PE) - z 

Here, z — exp(/3^i) is the fugacity with (3 = l/k B T and the chemical potential \i. The energy 
scale is chosen such that the minimum value of the d-dimensional external potential U(r) 
equals zero. Then the classical density of states for particles of mass m 

where the integration has to be taken over the classically accessible region, vanishes for negative 
energies. If p{E) — > for small energies, N reaches a finite maximum even for the maximum 
value of the fugacity z = 1. Additional particles arc then found in the ground state which is 
not included in the density of states (^). For a d-dimensional harmonic potential the density 
of states varies like E d ^ 1 and therefore Bose-Einstein condensation is only expected for d = 2 
and higher 0, |?| . 

This approach has been criticized by various authors. For the three-dimensional harmonic 
potential it was shown that corrections to the density of states (g) lead to corrections to the 
Bose temperature, which are of order iV" 1 / 3 [§. Other work § [fi| |ll|, || || has addressed 
the question whether there are effects due to the discreteness of the energy levels which is not 
included in (||). In this context, it was recently claimed that Bose-Einstein condensation can 
be found even in the case of a one-dimensional harmonic potential Q . 

The argument given by Ketterle and van Druten is as follows. The relative ground state 
occupation N /N is calculated by taking into account the discreteness of the energy levels. 
As an example, we show in fig. la results for particle numbers varying between N = 10 2 and 
N = 10 7 particles. The temperature is scaled with a critical temperature determined by 

»=¥vm ( 3) 



Tuo \ Too 

where hui is the energy level spacing. It can be seen that with increasing number of particles N, 
the crossover between finite and vanishing ground state occupation sharpens. This behaviour 
is similar to the three-dimensional case, where Bose-Einstein condensation occurs, and it is 
thus concluded that there is condensation even in one dimension. 

However, in a one-dimensional harmonic potential with fixed frequency u> the critical tem- 
perature T c diverges in the limit N — > oo. It is therefore not obvious how a phase transition in 
such a system may be defined. Furthermore, a rescaling with T c is then equivalent to widening 
the harmonic potential when increasing the particle number, i.e. lowering the oscillation 
frequency towards zero (cf. eq. (^)). 

We now address the question whether a macroscopic ground state population is a sufficient 
indicator to define Bose-Einstein condensation in this system. In this respect it is useful to 
consider another system with a constant density of states p(E). From eq. (^) we find for a 
two-dimensional box of length L 

P(E) = (4) 

27Tft 

Here a proper thermodynamic limit N, L — > oo, with constant density N/L 2 , can be performed 
and the notion of phase transition is well defined. While (2ir?i 2 / 'fcem)(JV '/ 'L 2 ) provides the tem- 
perature scaling, either N or L 2 determine how closely the thermodynamic limit is approached. 
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Fig. 1. - Relative ground state occupation No/N (a) for the one-dimensional harmonic po- 
tential and (b) for the two-dimensional box as a function of temperature T scaled with the 
critical temperature T c . The numerical results obtained by taking into account the discrete 
level structure correspond to N — 10 2 , 10 3 , 10 4 , 10 5 , 10 6 , and 10 7 particles from the upper right 
to the lower left curve. 



For the two-dimensional box it is generally accepted that, in contrast to the three-dimensional 
box, no phase transition occurs (see Q, where the authors derive this conclusion, although 
the specific heat for the two-dimensional box is not correctly reproduced |l4|). 

An analogy to the one-dimensional harmonic oscillator may be drawn by comparing its 
density of states p(E) = l/fiui to expression (|4|). The temperature scale is now given by NTiuj 
while the thermodynamic limit is reached either for N — > oo or hu — > 0. For fixed frequency to 
the limit N — > oo leads to a change in the temperature scale. A phase transition should there- 
fore become apparent when rescaling the temperature with N. For later purposes we note that 
the dimensionless density for the two-dimensional box N{\^/L) 2 , where At = {2-kTi 2 fi/m) 1 / 2 
is the thermal de Broglie wavelength, corresponds to Nfihtu for the one-dimensional harmonic 
oscillator. 

To check the analogy between the two systems, we have calculated the ground state pop- 
ulation for the two-dimensional box using discrete energy levels. Following the reasoning of 
Ketterle et al. we define a critical temperature via N = 2(L/At, c ) 2 ln(£/AT, c ) and obtain 
the results shown in fig. lb. Clearly, as in the case of a one-dimensional harmonic potential, 
the ground state population for large N shows a rather sharp transition to macroscopic values 
below the critical temperature. No qualitative difference in the behaviour of the two systems 
is visible. The definition of Bose-Einstein condensation via a macroscopic ground state pop- 
ulation enters therefore in conflict with predictions for the two-dimensional box. For a good 
criterion for condensation one instead has to resort to thermodynamic quantities. 

To this purpose we consider in particular the specific heat, which recently has become 
accessible experimentally |l5| 




(5) 



Here C is expressed through the continuous density of states p(E), which will prove to be 
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sufficient for our purposes. Evaluating (||) with constant density of states p = 1/Tud yields 

C 2 ^z™ f ln(l - z) dz 



Nk B Nfihuj ^ n 2 NHuj z d/3 

n—l 

and for the fugacity we obtain from 

2=1- cxp(-f3fujN). (7) 

This expression is correct for high temperatures but breaks down when z becomes of the order 
of its zero temperature value N/(N + 1), which is the case at temperatures of the order of 
the critical temperature defined by (||). At these temperatures z almost takes the value of 
one and the difference to the approximation (Q) is at most of order 1/N. Although (Q) could 
not be used to determine the ground state population at low temperatures, it allows us to 
obtain the correct result for the specific heat. Corrections of order 1/N to the fugacity are 
negligible in the first term on the right hand side of (||) whereas they are potentially dangerous 
in the logarithmic term. However in the temperature range of interest this term is strongly 
suppressed because dz/d/3 becomes exponentially small. 
With (0) the specific heat takes the form 

C = 2 ^ (1 - exp(-y))" y 
Nk B y 4r< n2 exp(y) - 1 

with y = phuN. The particle number N enters this result only through a linear rcscaling of the 
specific heat and the temperature. A similar result has been derived for a canonical ensemble 
[ fl2| . According to the above discussion the sum in (||) determines the low temperature 
behaviour of the specific heat which is then found to be 

Nk B 3 y 3N hui ' 1 ' 

In view of the analogy between the one-dimensional harmonic oscillator and the two-dimen- 
sional box the specific heat of the latter one then becomes 

C \ K 2 kv> ml? 



(10) 



In fig. 2, expression (|^) is depicted by the full line with its low temperature approximation 
([)[). The dashed line shows the result of a numerical calculation for 10 particles in a one- 
dimensional harmonic potential taking into account the discreteness of the spectrum. Even for 
such small N, (||) gives a rather good approximation. For 100 or more particles the difference 
between the result for the one-dimensional harmonic oscillator based on the discrete level 
structure on one hand and (|^) on the other hand becomes negligible. For the two-dimensional 
box the convergence of the calculation using discrete energies to the result for a constant 
density of states is somewhat slower due to the higher dimensionality of the system. 

Thus for large particle numbers the form of the specific heat becomes independent of N and 
neither singularities build up in the thermodynamic limit N — > oo nor a maximum appears 
at the critical temperature like the one found in the result sketched in S. Furthermore, 
the limiting result (^) is analytic for finite temperatures and we conclude that, as in the 
two-dimensional box, there is no phase transition in a one-dimensional harmonic potential. 

The difference between the behaviour of a non-interacting Bose gas in a one-dimensional and 
higher-dimensional harmonic potentials may also be understood in terms of discrete energy 
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Fig. 2. Fig. 3. 

Fig. 2. - Specific heat C of a Bose gas in a one-dimensional harmonic potential as a function 
of temperature T. The full line represents the result of a calculation with constant density 
of states with its low temperature approximation (dotted line), while the dashed line is for 
N — 10 particles taking into account the discrete level structure. For N = 10 eq. (||) yields 
kBT c /Nhu! = 0.454. Similar numerical results have been obtained in |jg| . 

Fig. 3. - Illustration of the reasoning leading to (|l5|). The lines correspond to the upper bound 
of the chemical potential ( |l2| ) and the high temperature approximation (|l4|), respectively. 



levels. An appropriate criterion for the possibility of Bose-Einstein condensation to occur is 
whether in the thermodynamic limit TV — > oo the chemical potential equals the ground state 
energy below the critical temperature or not. In the following we will show that this criterion 
is not met in a one-dimensional harmonic potential. We start out with the average number of 
particles in a grand canonical description 

i 00 (<*) 

N= - + V 9 - (11) 

exp(-/?/i)-l ^exp^ri-^-l' V ' 

where gffl is the degeneracy of the n-th eigenlevel of the d-dimensional harmonic potential. 



The first term on the right-hand-side of ( pL l[ ) represents the population of the ground state. 

One derives an upper bound of the chemical potential at arbitrary temperatures by consid- 
ering a situation where the excited states do not contribute to the number of particles 

a 1 knT 

i— < (12) 

Tito N Tilo 

Naturally this case coincides with the low temperature limit flhoj 3> 1. Chemical potential and 
temperature are taken with respect to the energy level spacing huj of the harmonic potential. 
This is the appropriate choice if u> is kept fixed in the thermodynamic limit. 

For high temperatures the ground state population becomes very small and therefore 
/3|/x| 3> 1. We may then approximate expression (|ll| ) for the average particle number by 

N = exp m e M -^n) = (13) 



Except for very low particle numbers the high temperature regime implies (3Tuo 1. The 
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chemical potential is then found as 

H _ In [N(j3huj) 



TlLO flflLO 



(14) 



and equals zero at hsT/tuo = N x / d . Since the ground state energy was taken to be zero, 
the high temperature result suggests that there might be Bose-Einstein condensation below a 
temperature TiuN 1 ' d /ks- A more detailed analysis leads to a prefactor of the order of one. 

To exclude the existence of a phase transition in a one-dimensional harmonic potential, it 
is sufficient to consider the upper bound (|l^) for the chemical potential at its intersection 
with the high temperature expression (0) . There the value of the chemical potential becomes 
(cf. fig. 3) 

JL = _jv(i-«0A». (15) 

For d = 1 this is of order one independently of the particle number N and one expects a 
smooth crossover at the intersection point. This behaviour forbids a phase transition in the 
one-dimensional case. On the other hand, for d > 2, the chemical potential at the intersection 
vanishes for N — ► oo opening the possibility of a sharp transition at a temperature of order 
N x / d , below which the chemical potential strictly equals the ground state energy. We thus 
recover the predictions of the theory based on the density of states (|^). 

In conclusion, by studying the specific heat and using an analogy with two-dimensional box 
systems we have shown that a macroscopic ground state population alone is not a sufficient 
indicator of a phase transition. In higher dimensional systems which show a condensation 
phenomenon a macroscopic ground state population is always accompanied by some signature 
in the specific heat or higher derivatives of the free energy. Their absence for the case of a 
one-dimensional harmonic potential leads us to conclude that no phase transition occurs in 
this system. The main interest of Bose-Einstein condensation lies in the particular coherence 
properties of condensates. Then the question arises to which extent the ground state popu- 
lation can give information about coherence. The results of this paper show that it might be 
necessary, maybe not sufficient, to carefully study thermodynamic quantities as well. 
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